Dynamically Driven Irreversibility in Many Body Quantum Systems: An 

experimental approach. 
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The NMR technique allows one to create a non-equilibrium local polarization and to detect its 
later evolution. By a change of the sign of the effective dipolar Hamiltonian, the apparently diffusive 
dynamics is reverted, generating a polarization echo (PE). This echo attenuates as a function of 
the time tu elapsed until the dynamics is reverted. We report a set of experiments where we slow 
down the dipolar dynamics showing that it controls this irreversibility. In cobaltocene, a cross over 
from Gaussian to exponential PE attenuation is found; in ferrocene the exponential regime is not 
reached. We conclude that many-body quantum systems have intrinsic instabilities, which amplify 
the environment fluctuations and hinder their backward evolution. 

PACS Number: 75.40. Gb, 76.60. Lz, 05.40. +j, 75.10.Jm. 



This letter reports the results of a set of NMR ex- 
periments designed to show the relationship between 
the complex nature of many-body dynamics and irre- 
versibility. Since Boltzmann's introduction jjj of the 
"Stosszahlansatz" the understanding of why an initial 
ordered state evolves irreversibly towards a higher en- 
tropy state has been a challenging issue 0. Loschmidt 
observed || that since microscopic laws are reversible, it 
should be possible to get a backward evolution recover- 
ing the initial state. However, this has remained in the 
realm of "gedanken" experiments, task of an imaginary 
Loschmidt's daemon, for almost a century. Recent theo- 
retical approaches to understand irreversibility point to 
the concept of local instabilities associated with classical 
chaotic systems (J] . If some initial state evolves accord- 
ing to a microscopic law for a time tp and then one makes 
the system evolve backwards, at time 2tp the state dif- 
fers from the original one by a "distance" proportional 
to exp [jtn] . The Lyapunov exponent 7 characterizes the 
chaotic motion responsible for the amplification J5| of any 
small error in the evolution. However, the quantum one- 
body counterpart Q], does not seem to present this in- 
stability leaving a conceptual gap in the formulation of a 
quantum dynamical description of irreversibility. Analyt- 
ical or numerical studies of many body quantum systems 
are lacking. 

Experimentally, a breakthrough into irreversibility was 
set up by E. Hahn with his NMR Spin Echo || which also 
applies to other spectroscopic techniques. It uses that the 
evolution of the macroscopic polarization of a system of 
N noninteracting spins, {-ft}, is described by the super- 
position of the one body Hamiltonians TL — ^ i TL{. An ex- 
ternal control parameter allows one to switch TLi — ► — TLi 
leading to a reversion of the dynamics. In this case, even- 
tual many body interactions limit the reversibility. Many 
body dynamics seemed much more difficult to invert until 
Rhim et al. achieved a Loschmidt's daemon for dipolar 
coupled spins. Since the sign and magnitude of this inter- 



action depends on the orientation of spins relative to their 
internuclear vector, the switch TL — * — TL is possible by 
changing the quantization axis. A further development 
is due to Zhang, Meier and Ernst § who were able to 
create and detect the local polarization of a spin Iy t repre- 
sented by the density matrix po = (5I + I\)/2 N ~ X . This 
allows one to deal with a well defined microstate |pj^0|- 
The recovered polarization after a forward and backward 
evolution is called Polarization Echo (PE). Later exper- 
iments pd| suggested that there are systems where the 
PE amplitude Mpe decays exponentially with time tp, 
while in others it decays with the much stronger Gaussian 
law. Figure la shows an idealized schematic view of the 
PE experiment. The initial state po evolves with TL + S, 
where £ describes small interactions not contained in the 
truncated dipolar Hamiltonian TL. The dashed curve in 
Fig. 1 is the local polarization detected at time t. It de- 
cays as 1 — d t 2 /2h 2 + 0(i 4 ), with d the dipolar second 
moment, evidencing an evolution according to the quan- 
tum many body dynamics ]l2fl . This polarization spreads 
beyond the neighborhood of site 1 in a typical time T m \, 
which is a few times h/d. If at a time tp we change 
the invertible part of the Hamiltonian TL — ► — TL, a PE 
builds up with a maximum {Mpe) at 2tp. The presence 
of uncontrolled processes prevents the perfect recovery 
of the initial state (Mpe < 1). Those "irreversible" in- 
teractions are simply the ones we do not know how to 
control. It is clear that if important interactions with 
a thermal bath were present they would show up as an 
exponential decay (Fig. 1) and S = A — ir. In this case 
the time w 7i/|r| is much shorter than r m b, and the 
irreversibility is directly controlled by the environment 
fluctuations. Our hypothesis is that, when r m b <C t^, the 
complex structure of the many body dynamics (N 3> 1) 
plays a crucial role in the attenuation of the PE. Then, 
intrinsic instabilities of the many body dynamics would 
lead to a stronger Gaussian decay (Fig. 1). This could 
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explain the amplification of fluctuations that, otherwise, 
could not produce an effective relaxation on the observed 
time scale. 

In order to test this hypothesis one must control the 
complexity of the state reached at time tp. This can be 
done in two ways: by changing the system (and hence the 
interaction network and relaxation mechanisms) or by re- 
ducing the effective dynamics. The latter provides a more 
controllable situation. Thus, in addition to selecting ap- 
propriate systems, we developed a pulse sequence which 
allows us to progressively reduce the complexity reached 
by the many body state, while the non-invertible inter- 
actions are kept constant. This is achieved by fitting n 
periods of alternating forward and backward dipolar evo- 
lutions in a fixed time tp. 

The sequence is schematized in Fig. 2. It uses the rare 
13 C spin (S'-spin) (1.1 % natural abundance) as a local 
probe that injects magnetization to one of the abundant 
1 H spins (/-spins) and later captures what is left. A 
(tt/2) x pulse on the /-spins creates a polarization that is 
transferred during tc to the S-spin when both are irra- 
diated at their respective resonant frequencies with field 
strengths fulfilling the Hartmann-Hahn condition. Af- 
ter the S-spin is polarized, it is kept spin-locked for a 
time ts while the coherence of the /-spins decays to zero. 
Then: A) A cross polarization (CP) pulse of duration 
td transfers the magnetization from an initially y polar- 
ized S'-spin to the y-axis of the /i-spin directly bonded to 
it |l3[ ], establishing the initial state po] B) The dipolar 
coupled /-spins evolve in presence of a strong spin-lock 
field that fixes the quantization axis and prevents the I-S 
coupling. If we neglect non-secular terms, the effective 
Hamiltonian 
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governs the forward evolution of po. The flip-flop terms 
1^ ' I// cause the spreading ( "diffusion" ) of the initial 
localized state; C) A (tt/2) pulse tilts the polarization 
to the laboratory frame where the /-spins evolve with 
— [2] H z while an S-irradiation prevents again the I-S 
coupling. A (tt/2)_ x pulse tilts the polarization back to 
the rotating xy plane. Taking the (it/ 2) pulses into ac- 
count, the effective Hamiltonian in this period is — [2] TL y . 
This change in the sign produces the backward evolution 
that builds up the PE. For a given time tR, we alternate 
n periods of forward and backward evolutions with du- 
rations of t\ = tp/n and r 2 = [i] T\ respectively. D) 
Another CP pulse of length is applied to transfer back 
the polarization to the x-axis of S. E) The S polarization 
is detected while the /-system is kept irradiated (high- 
resolution condition). The S-signal intensity is propor- 
tional to the remaining polarization in the /i-spin after 
the evolution periods. Since the S'-spin is not an ideal 



local probe (|Jl^], the proton spins have some forward 
evolution during the CP periods. To compensate for this 
undesired effect, it is necessary to introduce extra back- 
ward periods r' = [^] t^/2. In an ideal situation, the 
initial state should be recovered at the end of each cycle. 
The spreading of the local excitation due to the dipolar 
interaction and the number of dipolar spin correlations 
can be gradually reduced by increasing n. This sequence 
could be adapted to a broad range of polycrystals by 
complementing the magic angle spinning of the sample 
with rotor-synchronized pulses ]l4| . 

All the NMR measurements were made at room tem- 
perature in a Bruker MSL-300 spectrometer, equipped 
with a standard Bruker CP-MAS probe, operating at a 
13 C frequency of approximately 75.47 MHz. We ana- 
lyzed two different systems: a single crystal sample of 
ferrocene (C_5H 5 ) 2 Fe, and a polycrystalline sample 
of cobaltocene fll6fl , (C5H 5 ) 2 Co. The (C5H 5 ) 2 Fe experi- 
mental data were recorded for one of the two crystallo- 
graphically inequivalent sites, where the fivefold molec- 
ular symmetry axis is at 20° with respect to the exter- 
nal magnetic field. The (CsHs^Co data correspond to 
molecules with the symmetry axis perpendicular to the 
external field (they are frequency selected [jll| from the 
13 C spectrum). The relative phase and the rf amplitude 
of the four basic channels (X,—X,Y,—Y) are quite criti- 
cal and must be carefully adjusted $L7\ . 

Fig. 3 shows Mpe for the ferrocene sample as a func- 
tion of tR for n = 1, 2, 8, 16. We observe that its attenua- 
tion decreases with n and, in all the cases, the dynamical 
irreversibility manifests as a Gaussian decay. We obtain 
an excellent fitting (solid lines) with the characteristic 
time t™ 6 as the only free parameter (according to [ Q |, 
we required MpE(tR = 0) = 1 and an asymptotic value 
of 0). The inset shows that the dynamical time r mb has 
a linear dependence on n. The good signal to noise ratio 
in the single crystal for n = 1 allows us to ensure the 
Gaussian decay for over two orders of magnitude. This 
is a striking result since, from conventional wisdom, one 
could expect the attenuation of the PE to be directly 
caused by the non inverted interactions. These residual 
interactions, whatever they are, act for the same time tR 
independently of the value of n. Consequently, it is the 
reversible dipolar interaction the one which controls the 
decay of the PE in ferrocene. 

The situation of the cobaltocene sample is rather dif- 
ferent. Fig. 4 shows M PE (t R ) for n = 1,2,5,8,12. A 
crossover is evident between a dominant Gaussian at- 
tenuation to an exponential one. This can be explained 
as follows: While cobaltocene has the same interaction 
network as ferrocene, and therefore a similar dipolar dy- 
namics, it possesses a strong local source of relaxation 
(the paramagnetic Co atom) with a characteristic time 
t^. When n = 1, 2 the evolved state is complex enough, 
T mb < r 0' an d the dominant effect is the dynamical am- 
plification. However, when the dynamics is sufficiently 
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reduced (n = 8, 12, 16) one gets t™ 6 > and the ir- 
reversible mechanism becomes dominant. An asymp- 
totic pure exponential decay with ~ (970 ± 20) ps 
is reached upon further increase of n beyond 12. This 
paradigmatic example shows how the dynamical insta- 
bilities can be slowed down until the underlying mecha- 
nisms beyond our control are revealed. The crossover is 
not observed in ferrocene, since the time scale of the cor- 
responding exponential is beyond the experimental avail- 
ability. The solid lines in Fig. 4 represent fittings to 



the equation Mpe^r) = cxp 



Based on the results obtained for ferrocene we set r^ t 
an + b. In this way the whole set of data was fitted with 
only three free parameters. 

These experimental results may have important con- 
sequences in further developments of a theory of chaos 
in quantum many body systems. To see this, let pn+'E 
be the density matrix at time tp obtained from the 
evolution of po with a generalized master equation |is|] 
or the Keldysh equation (TJ)) for the field operators. 
Consider that further (backward) evolution is unitary, 
Uu = cxp(—ffl.tn/h). The normalized PE amplitude can 
be readily obtained: 



M PE {2t R ) = 2Tr (hU-Hpn+zUln 
= 2 N Tr(p H p n+ v)-l. 



(2) 



The last expression can be thought as the overlap be- 
tween two states evolved from the same initial state po 
with different Hamiltonians, 7i+£ and "H. From this 
point of view, the logarithm of Mpe has a suggestive 
resemblance to some entropy definitions [ ^0| . In semi- 
classical cases, an exponential decay is expected. In the 
non trivial case its rate 7 S> I/t^ depends only on 7i, 
and it is a manifestation of a strong sensitivity to slight 
changes in the dynamical parameters. This property is 
itself a signature of the presence of chaos |2lf| . However, 
our experimental result in a quantum many body system 
shows that the decay follows the much stronger Gaussian 
decay. This could be a consequence of the progressive 
availability of a huge Hilbert space made possible by the 
dynamics. 

The main conclusion of our experiments is that the 
complex dynamical structure presented by the many- 
body quantum systems has a very strong intrinsic in- 
stability towards irreversibility. Probably this instability 
was foreseen by the Boltzmann's deep physical insight. 
While we have a knob to slow down the many body dy- 
namics, whatever is left warrants an irreversible decay. 
Once we reach the time scale of a process we cannot con- 
trol it appears as an "ultimate" irreversible time. 

HMP acknowledges a fruitful discussion with Prof. E. 
L. Hahn at a Gordon Conference. HMP and PRL also 
benefited from an invitation to Zurich by Prof. R. R. 
Ernst at an early stage of this work and stimulating 



comments of Prof. A. Pines about our preliminary re- 
sults. This work was performed at LANAIS de RMN 
(UNC-CONICET), where the ferrocene single crystal was 
grown by Rodrigo A. Iglesias. Financial support was re- 
ceived from Fundacion Antorchas, CONICET, FoNCyT, 
CONICOR and SeCyT-UNC. The authors are affiliated 
to CONICET. 



[8. 

[9 
[10 

[11 

[12; 

[13; 

[14 
[15; 
[16 
[17 

[is; 

[19 

[20 



L. Boltzmann, Ann. Phys. (Leipzig) 57, 773 (1896); 
translated and reprint in S. G. Brush, Kinetic Theory 
2, Pergamon, Elmsford (1966). 

J. L. Lebowitz, Physics Today, 46 (9) 32 (September 
1993); J. J. Halliwell, J. Perez-Mercader and W. H. 
Zurek, eds.. Physical Origins of Time Asymmetry, Cam- 
bridge Univ. Press. (Cambridge 1994). 
J. Loschmidt: Sitzungsber. Kais. Akad. Wiss. Wien 
Math. Naturwiss. Klasse 73, 128-142 (1876). 
F. M. Izrailev, Phys. Rep. 196, 299 (1990); G. Casati, B. 
V. Chirikov, D. L. Shepelyansky and I. Guarneri, Phys. 
Rep. 154, 87 (1987). 

A. Bellemans and J. Orban, Phys. Lett. 24A, 620 (1967). 

E. L. Hahn, Phys. Rev. 80, 580 (1950); R. G. Brewer and 

E. L. Hahn, Scientific American, Dec. 1984 p50. 

W. K. Rhim, A. Pines and J. S. Waugh, Phys. Rev. Lett., 

25 , 218 (1970); Phys. Rev. B 3, 684 (1971); W. K. Rhim 

and H. Kessemeier, Phys. Rev. B 3, 3655 (1971). 

S. Zhang, B. H. Meier, and R. R. Ernst, Phys. Rev. Lett. 

69, 2149 (1992); Solid State Nucl. Magn. Reson. 3, 121 

(1994). 

H. M. Pastawski, PR. Levstein, and G. Usaj, Phys. Rev. 
Lett. 75, 4310 (1995). 

H. M. Pastawski, G. Usaj and PR. Levstein, Chem. 
Phys. Lett. 261, 329 (1996). 

P. R. Levstein, G. Usaj and H. M. Pastawski, J. Chem. 
Phys. (A7.06.021 at press). 

H. M. Pastawski and G. Usaj, Phys. Rev. B (BBJ643 at 
press) . 

L. Muller, A. Kumar, T. Baumann, and R. R. Ernst, 
Phys. Rev. Lett. 32 1402 (1974). 

M. Tomaselli, S. Hediger, D. Suter and R. R. Ernst, J. 
Chem. Phys. 105, 10672 (1996). 

P. Seiler, and J. D. Dunitz, Acta Crystallog. B 35 1068 
(1979). 

M. Y. Antipin, R. Boese, N. Auyart and G. S. Schmidt, 
Structural Chem. 4, 91 (1993). 

D. P. Burum, M. Linder and R. R. Ernst, J. Magn. Res. 
43, 463 (1981). 

R. R. Ernst, G. Bodenhausen and A. Wokaun: Principles 
of Nuclear Magnetic Resonance in One and Two Dimen- 
sions (Clarendon, Oxford, 1987). 

L. V. Keldysh, Zh. Eksp. Teor. Fiz 47, 1515 (1964) [Sov. 
Phys. JETP 20, 1018 (1965)] ;H. M. Pastawski, 46 4053 
(1992). 

A. K. Pattanyak and P. Brumer, Phys. Rev. Lett. 79, 
4131 (1997). 



3 



[21] F. Haakc, Quantum Signature of Chaos Springer- Verlag 
(1992). 

.Figure 1: Left panel: Idealized representation of a 
Polarization Echo experiment. The local polarization de- 
cay (dotted line) is described by the Hamiltonian H + S, 
where X is small enough for TL to determine the spreading 
time T m b- At time tft the dynamics is reversed. Further 
evolution of the polarization is shown by the solid line. It 
develops a PE of amplitude Mpe- Right panel: Possi- 
ble dependences of Mp e as a function of tp . The dotted 
line represents a system controlled by a strong relaxation 
rate (t^ <C r m t) while the solid line shows the opposite 
case. The time scale r in the plot corresponds to and 
r m fc respectively. 

Figure 2: Pulse sequence to control the dipolar spin 
dynamics after creating a local polarization (see text). 
The solid and dotted lines show the main and secondary 
amplitude polarization pathways. 

Figure 3: Attenuation of the polarization echo in fer- 
rocene as a function of tp. The data were recorded using 
the pulse sequence schematized in Fig. 2 with: tc = 2ms, 
ts = 1ms, td = 53/is, and n = 1,2,8,16. The attenu- 
ation is slowed down as the complexity reached by the 
system is reduced. The lines represent Gaussian fittings 
with the characteristic time T^ b as the only free param- 
eter. Inset: Characteristic time dependence on n. 

Figure 4: Same as in Fig. 3 but for cobaltocene sam- 
ple. We used tc = 85fis, ts = 150yus, td = 85fis, and 
n = 1,2,5,8,12. There is a clear crossover between a 
Gaussian and an exponential decay. This crossover re- 
sults from a change of the relative importance between 
r^b an d i>. The solid lines represent fittings of the whole 
set of data with only three free parameters (see text). 
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